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Abstract. We investigate the role which clouds could play in 
resolving the Faint Young Sun Paradox (FYSP). Lower solar 
luminosity in the past means that less energy was absorbed 
on Earth (a forcing of — 50 Wm~ 2 during the late Archean), 
but geological evidence points to the Earth being at least as 
warm as it is today, with only very occasional glaciations. 
We perform radiative calculations on a single global mean 
atmospheric column. We select a nominal set of three lay- 
ered, randomly overlapping clouds, which are both consis- 
tent with observed cloud climatologies and reproduce the ob- 
served global mean energy budget of Earth. By varying the 
fraction, thickness, height and particle size of these clouds 
we conduct a wide exploration of how changed clouds could 
affect climate, thus constraining how clouds could contribute 
to resolving the FYSP. Low clouds reflect sunlight but have 
little greenhouse effect. Removing them entirely gives a forc- 
ing of +25 W m~ 2 whilst more modest reduction in their effi- 
cacy gives a forcing of +10 to +15 Wm~ 2 . For high clouds, 
the greenhouse effect dominates. It is possible to generate 
+50Wm~ 2 forcing from enhancing these, but this requires 
making them 3.5 times thicker and 14 K colder than the stan- 
dard high cloud in our nominal set and expanding their cover- 
age to 100% of the sky. Such changes are not credible. More 
plausible changes would generate no more that +15 Wm -2 
forcing. Thus neither fewer low clouds nor more high clouds 
can provide enough forcing to resolve the FYSP. Decreased 
surface albedo can contribute no more than +5 Wm~ 2 forc- 
ing. Some models which have been applied to the FYSP do 
not include clouds at all. These overestimate the forcing due 
to increased C0 2 by 20 to 25% when pC0 2 is 0.01 to 0. 1 bar. 



1 Introduction 

Earth received considerably less energy from the Sun early 
in Earth's history than today; ca. 2.5 Ga (billion years be- 



fore present) the sun was only 80% as bright as today. Yet 
the geological evidence suggests generally warm conditions 
with only occasional glaciation. This apparent contradiction 



is known as the Faint Young Sun Paradox (FYSP, Ringwood 



|1961| Sagan and Mullen 1972| l. A warm or temperate cli- 
mate under a faint sun implies that Earth had either a stronger 
greenhouse effect or a lower planetary albedo in the past, or 
both. In this study, we focus on the role of clouds in the 
FYSP. We both examine how their representation in models 
affects calculations of changes in the greenhouse effect, and 
constrain the direct contribution that changing clouds could 
make to resolving the FYSP. 

Clouds have two contrasting radiative effects. In the spec- 
tral region of solar radiation (shortwave hereafter) clouds 
are highly reflective. Hence clouds contribute a large part 
of Earth's planetary albedo (specifically the Bond albedo, 
which refers to the fraction of incident sunlight of all wave- 
lengths reflected by the planet). In the spectral region of 
terrestrial thermal radiation (longwave hereafter), clouds are 
a strong radiative absorber, contributing significantly to the 
greenhouse effect. Cloud absorption is largely independent 
of wavelength (they approximate to "grey" absorbers), in 
contrast to gaseous absorbers which absorb only in certain 
spectral regions corresponding to the vibration-rotation lines 
of the molecules. 

Despite the obvious, first-order, importance of clouds in 
climate, it has become conventional to omit them in mod- 
els of early Earth climate and use instead an artificially high 



surface albedo. As described by Kasting et al. ( 1984): 
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Clouds are not included explicitly in the model; 
however, their effect on the radiation balance is 
accounted for by adjusting the effective albedo to 
yield a mean surface temperature of 288 K for the 
present Earth. The albedo is then held fixed for 
all calculations at reduced solar fluxes. ...we feel 
that the assumption of constant albedo is as good 
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as can be done, given the large uncertainties in the 
effect of cloud and ice albedo feedbacks. 

In effect, the surface is whitewashed in lieu of putting clouds 
in the atmosphere. This assumption has been used exten- 



sively in the models from Jim Kasting's group ( 


lasting et al. 


1984 Kasting and Ackerman 


1986; Kasting 


1987||1988 


Kasting et al.| 


1993[|Pavlov et al. 


2000, 2003 


Kasting and 


Howard 2006; Haqq-Misra et al.| |2008 1, which, together 



with parametrisations and results based on these models (for 



example, Kasting et al.||1988 ; Caldeira and Kasting , 1992a b; 


Kasting et al. 


2001; Kastingr|2005t |Tajika| |2003T|Bendt- 


sen and Bjerrum 




2000; von Bloh etaL] |2003a|b| [Lenton and von Bloh |200T] 



Bergman et al.| |2004[ ) have dominated early Earth palaeocli- 



mate and other long term climate change research for the last 
two and a half decades. The validity of this method has not 
previously been tested. 

Whilst Kasting's approach is that changes to clouds are so 
difficult to constrain that one cannot justifiably invoke them 
to resolve the FYSP, others are more bold. Some recent pa- 
pers have proposed cloud-based resolutions to the FYSP. 
Rondanelli and Lindzen| ( |2010[ ) focus on increasing the 



warming effect of high clouds, finding that a total covering 
of high clouds which have been optimised for their warm- 
ing effect could give a late Archean global mean temper- 
ature at freezing without increasing greenhouse gases. To 
justify such extensive clouds, they invoke the "iris" hypoth- 
esis ( |Lindzen et al. 2001 1 which postulates that cirrus cov- 
erage should increase if surface temperatures decrease (this 



hypothesis has received much criticism, e.g. Hartmann and 
|Michelsen| [20021 |Chambers et al.||2002) . 

|Rosing et al.| ( |2010| l focus on decreasing the reflectivity of 
low level clouds so that the Earth absorbs more solar radia- 
tion. To justify this, they suggest that there was no emission 
of the important biogenic cloud condensation nuclei (CCN) 
precursor dimethyl sulphide (DMS) during the Archean and, 
consequently, clouds were thinner and had larger particle 
sizes. 

We note that both |Rondanelli and Lindzen|(2010| ) and |Ros-| 
|ing et al7] ( |2010| l predict early Earth temperatures substan- 
tially below today's, which we do not consider a satisfactory 
resolution of the FYSP. 



Shaviv ( 2003 i and Svensmark ( 2007 1 proposes less low 



level cloud on early Earth due to fewer galactic cosmic 
rays being incident on the lower troposphere. The under- 
lying hypothesis is of a correlation between galactic cos- 
mic ray incidence and stratus amount, through CCN cre- 



ation due to tropospheric ionization (Svensmark and Friis- 
Christensen 1997; Svensmark 2007]). This has received ex- 
tensive study in relation to contemporary climate change and 



has been refute d (e.g.|Sun and Bradley||2002||Lockwood and| 
Frohhch|[2007l|Kristjansson et al.||2008[|Bailer-Jones| |2009[ 
Calogovic et al.||2010||Kulmala et al.||2010| and references 



In this study, we comprehensively asses how the radia- 
tive properties of clouds, and changes to these, can affect the 
FYSP. First, we explicitly evaluate how accurate cloud-free 
calculations of changes in the greenhouse effect are with re- 
spect to atmospheres with clouds included. We do this by 
considering a very wide range of cloud properties within 
a single global mean atmospheric column, finding a case 
study which matches Earth's energy budget, then comparing 
the effect of more greenhouse gas in this column to a cloud- 
free calculation. Second, we conduct a very wide exploration 
of how changing clouds could directly influence climate. We 
vary fraction, thickness, height and particle size of the clouds 
and vary surface albedo. We do not advocate any particular 
set of changes to clouds. Rather, constraints on what direct 
contribution clouds can make towards resolving the FYSP 
emerges from our wide exploration of the phase space. 

The heyday of clouds in 1-D models was the 1970s and 
1980s. Improvements in radiative transfer codes and compu- 
tational power over the last 30 years allow us to contribute 
new insight to the problem, in particular by widely exploring 
phase-space. Nonetheless, these classic papers retain their 
relevence and are instructive as to how one might treat clouds 
in such simple models (e.g. |Schneider| 1972 Reck 



|Wang and Stone] [1980} [Stephens and Websterj 1981 



lock||1982 i. With specific relevance to the FYSP, |Kiehl and 



1979 



Char- 



|Dickin son ( 1987 ) included clouds in their model of methane 
and carbon dioxide warming on early Earth, and calculated 
radiative forcings from some changed cloud cases (our re- 
sults here agree with this older work). |Rossow et al. ( 1982| l 
considered cloud feedbacks for early Earth. 

Regarding whether a cloud-free model will correctly cal- 
culate the increased greenhouse effect with increases gaseous 
absorbers, we hypothesise that it will lead to an overestima- 
tion in the efficacy of enhanced greenhouse gases. In the 
absence of clouds, the broadest range of absorption is due to 
water vapour. However, whilst water vapour absorbs strongly 
at shorter and longer wavelengths, it absorbs weakly between 
8 and 15 um. This region of weak absorption is known as the 
water vapour window. It is coincident with the Wein peak of 
Earth's surface thermal emission at 10 um. Thus the water 
vapour window permits a great deal of surface radiation to 
escape to space unhindered. Other greenhouse gases - and 
clouds - do absorb here, so are especially important to the 
greenhouse effect. With clouds absorbing some fraction of 
the radiation at all wavelengths, the increase in absorption 
with increased greenhouse gas concentration would be less 
than if clouds were absent. Therefore, we think that a cloud- 
free model would overestimate increased gaseous absorp- 
tion with increased greenhouse gas abundance and underes- 
timate the greenhouse gas concentrations required to keep 
early Earth warm. 

Comparison of cloudy and cloud-free radiative forcings in 
the context of anthropogenic climate change (Pinnock et al. 



therein). 



1995] |Myhre and Stordalj [1997] |Jain et al.| [2000) supports 



our hypothesis. For CO2, a clear-sky calculation overesti- 
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mates the radiative forcing by 14%. For more exotic green- 
house gases, which are optically thin at standard conditions 
(CFCs, CCs, HCFCs, HFCs, PFCs, bromocarbons, iodocar- 
bons), clear-sky calculations overestimate radiative forcing 
by 26-35%. CH4 and N2O are intermediate; their clear-sky 
radiative forcings are overestimated by 29% and 25%, re- 
spectively ( |Jain et aL] 2000). 



Table 1. GAM profile at levels (layer boundaries). Note that the 
tropopause is at lOOhPa. 



A roadmap of our paper is as follows. In Sect. [2] we de- 
scribe our general methods, verification of the radiative trans- 
fer scheme and the atmospheric profile we use. In Sect. [3] 
we deal specifically with the development of a case study of 
three cloud layers representing the present climate and the 
model sensitivity to this. In Sect. [4] we compare cloudy and 
cloud-free calculations of the forcing from increased green- 
house gas concentration. In Sect. [5] we explore what direct 
forcing clouds could impart, and in section [6] we evaluate 
the aforementioned cloud-based hypotheses for resolving the 
FYSP. 



2 Methods 
2.1 Overview 

Using a freely available radiative transfer code, we develop 
a set of cloud profiles for single-column models which is in 
agreement both with cloud climatology and the global mean 
energy budget. This serves as the basis for comparison of 
radiative forcing with a clear sky model and with changed 
cloud properties. 



2.2 Radiative forcing 



Pressure 


A Iti til r\f* 


Teniperciture 


^Viiter vapour 


Ozone 




V Km / 






(ppmv) 


10 


64.739 


230.00 


0.0036 


1.080 


20 


59.912 


245.61 


0.0036 


1.384 


30 


56.951 


252.88 


0.0036 


1.626 


50 


53.114 


260.00 


0.0036 


1.974 


100 


47.763 


266.29 


0.0035 


2.600 


200 


42.393 


260.33 


0.0033 


5.484 


300 


39.339 


254.31 


0.0032 


6.810 


500 


35.612 


243.24 


0.0032 


7.242 


1000 


30.842 


228.11 


0.0031 


7.490 


2000 


26.290 


222.10 


0.0030 


6.169 


3000 


23.671 


218.71 


0.0029 


4.780 


5000 


20.445 


212.59 


0.0026 


2.250 


10000 


16.204 


206.89 


0.0023 


0.516 


15000 


13.727 


211.83 


0.0048 


0.344 


20000 


11.914 


219.01 


0.0153 


0.160 


25000 


10.461 


225.87 


0.0456 


0.122 


30000 


9.237 


233.27 


0.1852 


0.089 


35000 


8.168 


240.52 


0.3751 


0.070 


40000 


7.215 


247.19 


0.6046 


0.058 


45000 


6.352 


253.27 


0.8866 


0.051 


50000 


5.562 


258.62 


1.2365 


0.047 


55000 


4.834 


263.15 


1.6525 


0.045 


60000 


4.159 


267.14 


2.1423 


0.045 


65000 


3.529 


270.73 


2.7049 


0.044 


70000 


2.938 


274.00 


3.3366 


0.042 


75000 


2.381 


277.05 


4.1602 


0.039 


80000 


1.855 


279.84 


5.2152 


0.035 


85000 


1.356 


282.28 


6.3997 


0.033 


90000 


0.882 


284.08 


7.8771 


0.032 


95000 


0.431 


285.85 


9.5702 


0.031 


100000 


0.000 


289.00 


11.1811 


0.031 



In work on contemporary climatic change, extensive use is 
made of radiative forcing to compare the efficacy of green- 
house gases (e.g. Forster et al. 20071. This is defined as 



the change in the net flux at the tropopause with a change in 
greenhouse gas concentration, calculated either on a single 
fixed temperature-pressure profile or a set of fixed profiles 
and in the absence of climate feedbacks. Surface tempera- 
ture change is directly proportional to radiative forcing, with 
a radiative forcing of approximately 5 W vaT 2 being required 
to cause a surface temperature change of 1 K (see Fig. 7 of 
2009b). Note that the tropopause must be 
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defined as the level at which radiative heating becomes the 
dominant diabatic heating term (Forster et al. 1997| ), i.e. the 
lowest level at which the atmosphere is in radiative equilib- 
rium. 

We base all our analyses on radiative forcings here. As 
we make millions of radiative transfer code evaluations, sav- 
ing in computational cost from comparing radiative forcings 
rather than running a radiative-convective climate model is 
significant, and facilitates the wide range of comparisons pre- 
sented. 



2.3 Global Annual Mean atmosphere 

We perform all our radiative transfer calculations on a sin- 
gle Global Annual Mean (GAM) atmospheric profile (Ta- 
ble [TJ. This is based on the GAM profile of Christidis et al. 
(1997) with some additional high altitude data from Jain 
let al.|(| 2000] >. Surface albedo is set as 0.125 (Trenberth et al. 



[2009 



). For standard conditions we use year 2000 gas compo- 
sitions: 369 ppmv C0 2 , 1760 ppbv CH 4 and, 316ppbv N 2 0. 
We use present day oxygen and ozone compositions through- 
out the work. Whilst comparisons without these might be 
interesting, they would necessitate using a different temper- 
ature profile in order to be self consistent. This would signif- 
icantly complicate our methods, so no such calculations are 
performed. For solar calculations, we use the present solar 
flux and a zenith angle of 60° . 

Calculating radiative forcings on a single profile does in- 
troduce some error relative to using a set of profiles for var- 



ious latitudes (Myhre and Stordal 1997 Freckleton et al. 



1998 1 Jain et al. 2000 ). However, as this is a methodological 



4 



C. Goldblatt and K. J. Zahnle: Clouds and the Faint Young Sun Paradox 



paper concerning single column radiative-convective models, 
it is the appropriate approach to take here. 

2.4 Radiative transfer code and verification 

We use the Atmosphere Environment Research (AER) Rapid 
Radiative Transfer Model (RRTM, |Mlawer et IE] \TW7\ 
Clough et al. 2005), longwave version 3.0 and shortwave 



version 2.5, which are available from http://rtweb.aer.com 
(despite different version numbers, these were both the most 
recent versions at the time of the research). RRTM has been 
parameterised for pressures between 0.01 and 1050 hPa and 
for temperatures deviating no more than 30 K from the stan- 
dard mid-latitude summer (MLS) profile. We have verified 
that the GAM profile we use is within this region of pressure- 
temperature space. The cloud parameterisations in RRTM 
which we select follow Hu and Stamnes| (1993 ) for water 
clouds and |Fu et al.| ( |1998] l for ice clouds. 

RRTM has been designed primarily for contemporary at- 
mospheric composition. Our intended use is for different 
atmospheric composition (higher greenhouse gas concentra- 
tions), so it is necessary for us to independently test the 
performance of the model under these conditions ( |Collins| 
|et~aTl [20061 |Goldblatt et al.| |2009b| >. Following the ap- 
proach of Goldblatt et al. (2009b) we directly compare long- 
wave clear sky radiative forcings from RRTM to the AER 



Line-by-Line Radiative Transfer Model (LBLRTM, Clough 



|et al.||2003] ). These runs are done on a standard Mid-Latitude 
Summer (MLS) profile ( |McClatchey et al.| |1971[ |Anderson| 



Our focus is on comparison of cloudy to cloud-free pro- 
files within RRTM, so we do not need high accuracy cal- 
culations of early Earth radiative forcings. We can there- 
fore use rather relaxed and qualitative thresholds for ac- 
ceptable model performance relative to LBLRTM: we re- 
quire continuous and monotonic response to changing green- 
house gas concentration (no saturation), the forcing should 
be smooth and monotonic and divergence from the LBLRTM 
flux gradient should be limited. For CO2, RRTM forcing 
is not smooth or monotonic below pCO2=10~ 4 bar so this 
region is excluded (see Fig. [TJ. CO2 concentrations up to 
pCO2=10~ 1 bar are used, though there is some underestima- 
tion of radiative forcing by RRTM above pCO2=10 -2 bar. 
Also, collision-induced absorption (absorption due to for- 
bidden transitions) becomes important at J3CO2 ~ 0. 1 bar (J. 
Kasting, private communication) but coefficients for these 
are not included in the HITRAN database on which both 
RRTM and LBLRTM absorption coefficients are based. 
Therefore, it is emphasised that the radiative forcings pre- 
sented here for high CO2 will be underestimates, but valid 
for intra-comparison. 

The comparison of RRTM to LBLRTM (Fig.[T) is only for 
the purpose of validating clear sky radiative forcing in the 
context of this methodological study. We have undoubtedly 
used RRTM outside its design range. This is not intended as 
an assessment of its use for the contemporary atmosphere or 
for anthropogenic climate change. 



etaEl |1986) to take advantage of the large number o f com- 3 cloud representation and model tuning 



putationally expensive LBLRTM runs performed by Gold- 



|blatt et aT| ( |2009r5] >. Performance of the codes is evaluated 
at three levels: the top of the atmosphere (TOA), the MLS 
tropopause at 200 hPa and the surface. Upward and down- 
ward fluxes are considered separately. The surface is taken 
to be a black body, so the upward flux depends only on tem- 
perature (F^ w surf = aT 4 ). The downward longwave flux at 
the TOA is zero. Neither vary with greenhouse gas concen- 
trations, so changes in the net flux at these levels depends on 
one radiation stream only. At the tropopause the net flux is 
the sum of the two streams. It is defined positive downwards, 



Pi — F^ — F^ 



(1) 



In addition to the radiative flux, we show (Fig.[T]i the forcing 



F 



lw 



_Fi w — F\ 



lw.std ) 



(2) 



where -Fiw.std is the flux at preindustrial conditions and the 
flux gradient (change of flux with changing gas concentra- 
tion) 



dF 
dX 



AF 
AX 



Fi+i — Fi 



(3) 



3.1 Practical problems and observational guidance 

Generation of an appropriate cloud climatology for this work 
is not straightforward. Two fundamental problems are short- 
comings in available cloud climatologies and averaging to 
a single profile. Concerning climatologies, the problem is 
one of observations: surface observers will see the lowest 
level cloud only, satellites will see the highest level of cloud 
only. Radiosondes are cloud penetrating and cloud proper- 
ties may be inferred from measured relative humidity, but the 
spatial and temporal coverage of radiosonde stations is lim- 
ited. See [Wang etaE] ( |2000| l and |Rossow et al.| ( |2005] l for ex- 
tensive discussion of what progress can be made. Similarly, 
radar can profile clouds, but such observations are sparse. 
Concerning averaging, the dependence of the global energy 
budget on cloud properties is expected to be non-linear: one 
should not expect that a linear average of global cloud phys- 
ical properties would translate into a set of clouds whose ra- 
diative properties would give energy balance in a single co- 
lumn. Nonetheless, available temporally and spatially aver- 
aged data for cloud properties can guide how clouds should 
be represented in the model. 



where Ff is the flux at gas concentration X^ (Goldblatt et al. 
|2009b) . 



Rossow et al. ( 2005 1 deduce zonally-averaged cloud frac- 



tion profiles using a combination of International Satellite 
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Fig. 1. RRTM performance for CO2 for each flux and each level. Colours (online only) and markers are: black + for LBLRTM, magenta 
x for RRTM. Shaded areas are range from Quaternary minimum (180ppmv) to SRES maximum (1248 ppmv) concentration. Grey lines in 
these areas are solid for pre-industrial (287 ppmv) and dashed for year 2000 (369 ppmv) concentrations. 
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Fig. 2. Average cloud fraction with altitude following Rossow et al. 
( |2005| and W. Rossow, personal communication, 2009), for January 
and July, land and ocean. White areas are where there is either no 
land (the Southern and Arctic Oceans) or no ocean (Antarctica). 



Cloud Climatology Program (ISCCP) and radiosonde data 
(Fig. [2]). The existence of three distinct cloud layers and the 
pressure levels of these are immediately apparent when av- 
eraging the data meridionally (Fig. [3}. Following |Rossow| 
|and Schiffer| ( |1999| l, we divide the clouds into three groups, 
with divisions at 450 and 700 hPa. "Low" cloud coresponds 
to cumulus, stratocumulus and stratus clouds. "Mid" level 
clouds correspond to altocumulus, altostratus and nimbo- 
stratus. "High" clouds correspond to cirrus and cirrostra- 
tus. Absolute cloud fractions cannot be extracted directly 
from these data as information on how the clouds overlap 
is lost in temporal and spatial averaging. A simple ap- 
proach to give indicative values is to assume either maxi- 
mum or random overlap within each group (high, middle and 
low), then to scale these cloud amounts by a constant such 
that randomly overlapping the three groups gives the IPCC 



mean global cloud fraction of 67.6% (Rossow and Schif- 
fer 1999[ ). Maximum and random overlap within groups 
give cloud fractions [/high, /mid, /low] = [0. 24, 0.25, 0.43] and 
[/high, /mid , /low] = [0.25,0.29,0.39], respectively. 

Averaged cloud optical thickness or water paths are more 
difficult to constrain, as they are not directly available 
from the |Rossow et al.| ( |2005| l data set (W. Rossow, per- 
sonal communication, 2009). We proceed with ISCCP data 
only. Rossow and Schiffer| ( |1999 1 report water paths of 
[W h igh,Wmid,Wi O w]=[23,60,51]gm- 2 . ISCCP data are 
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Fig. 3. Average cloud fraction with altitude following Rossow et al. 
( |2005| and W. Rossow, personal communication, 2009). (a) North- 
ern Hemisphere, (b) Southern Hemisphere and (c) global for Jan- 
uary (green), July (purple) and mean (black). Grey horizontal lines 
separate high, mid- and low-level clouds. 



from downward looking satellite data only and overlap is not 
accounted for. Whilst the low cloud value will indicate low 
clouds only, high and mid level cloud values may include 
opacity contributions from the lower clouds which they ob- 
scure (see Fig. |2j. Hence these water paths are indicative 
only. 

Using fine resolution spatially and temporally resolved 
data would help resolve these issues. However, to do so 
would be beyond the scope of this work and, we feel, it is 
beyond what is necessary to address the first-order questions 
which are the subject of this paper. 

3.2 Development of cloud profiles 

We need to develop a set of cloud profiles which appropri- 
ately represents Earth's cloud and energy budget climatolo- 
gies. By necessity, we shall need to simplify cloud proper- 
ties, tune our model clouds and consider sensitivity of the 
model energy budget to these clouds. 



Table 2. Parameter values used in the large cloud tuning ensemble. 
Optical depth depends logarithmically on water path. Water path 
and cloud fraction are varied independently for each layer. Water 
paths range from optically thin to optically thick clouds jCurry arid| 
|Webster||1999} with 10 values. Cloud fractions range from 5% to 
100% coverage with 20 cases. Effective radius is for water clouds 
(low and mid level) and generalised effective size is for ice clouds 
(high). There there are 10 3 x 20 3 =8x 10 6 cases in total. 



Fixed properties 


High 


Mid 


Low 


Cloud top (hPa) 


300 


550 


750 


Cloud base (hPa) 


350 


650 


900 


Liquid or ice 


Ice 


Liquid 


Liquid 


Effective radius (pm) 




11 


11 


Generalised effective size (urn) 


75 






Variable properties 




All layers 




Water path (gm~ 2 ) 


[lO 04 


,10 06 ,10 a8 , 


..,10 2 ' 2 ] 


Cloud fraction 


[0.05,0.10, 0.15, . 


., 1.00] 



Even with the assumption that each cloud is homoge- 
neous, each of our three cloud layers is represented by 
a cloud base and top, water path, liquid:ice ratio, and 
effective particle sizes for liquid and ice particles, giv- 
ing 6 degrees of freedom for each cloud. With three lay- 
ers, there are eight permutations for overlap, contributing 
another 7 degrees of freedom for the fractional coverage. 
A total of 25 degrees of freedom is clearly impossible to 
explore fully. As a necessary simplification, we fix the 
cloud base and top, take clouds to be either liquid (low 
and mid clouds) or ice (high clouds) and fix the particle 



size (following Rossow and Schiffer 1999 1. We assume that 
cloud layers are randomly overlapped, so each cloud layer 
can be represented by a single fraction from which the over- 
lap is calculated (many GCMs use a "maximum-random" 
overlap method where cloud fractions in adjacent layers are 
correlated; this is not relevant here as our discrete cloud lay- 



ers are separated by intervals of clear sky, e.g. see Hogan and 
mingworth||2000) . 



Random overlap is easiest to explain for the case of two 
cloud levels (A and B), with cloud fractions a and b. Fraction 
ab of the sky would have both cloud layers, fraction o(l— b) 
would only have level A clouds and fraction (1— a)b would 
only have level B clouds, and fraction (1— a) (1— b) would be 
cloud free. With three cloud layers, we have eight columns. 
Each column is evaluated separately in both longwave and 
shortwave spectral regions and the final single column is 
found as a weighted sum of these devaluations. Different 
cloud fractions can be accounted for in this summation, re- 
ducing the number of RRTM evaluations needed. 

For each cloud layer, cloud fraction and water path are 
varied widely whilst the other four parameters are fixed (Ta- 
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ble|2|i. In all of the resulting cloud cases, we run the radia- 
tive transfer code for both standard and elevated CO2 levels 
(369 ppmv and 50 000 ppmv), giving 16 million runs in total. 

We refer to model runs in which we include clouds in this 
way as including "real clouds". This is meant in the sense 
that clouds are included in the radiative transfer code in a 
detailed and physically based manner. This is by contrast to 
previous models (e.g. Kasting et al.| 1984) >, where clouds are 



represented non-physically by changing surface albedo. 
3.3 Sensitivity experiment 

For each cloud case that we have defined (the large ensemble, 
Table[2]), we calculate the radiative forcing at the tropopause 
(F trop ), to which change in surface temperature is propor- 
tional. Radiative forcing is the change in net flux (here with 
increasing CO2): 



Ftrop — -F]trop,highC0 2 ] — -F[trop,stdC0 2 



(4) 



where F in each case is the net flux as a sum of longwave 
and shortwave spectral regions and upward and downward 
streams of radiation: 



SW/ 



(5) 



We consider two subsets of the large ensemble: 



1 . Cloud sets which give energy balance at the TOA. This 
is the most basic constraint on a possible climate. With 
|-Ftoa|<5 Wm~ 2 , a subset of 1.0 million cases re- 
mains. A relatively large |i 7 TOA,stdC0 2 1 i s allowed as 
variations in the water path are coarse, but it is corrected 
for by calculating radiative forcings so cannot bias the 
outcome. 



Cloud sets which give energy balance at the 
TOA and are close to observed longwave and 
shortw ave fluxes at the TOA ( |Trenberth et al."] 
2009). Constraints are |F TO A,stdC0 2 |<5 Wm~ 2 , 
95<F tTOA,sw,stdco 2 <115Wm- 2 
227<F tTOA.LW.stdco, <247Wm- 2 ). This 



I TOA,LW,stdCO 

a subset of 36 985 cases 



and 

gives 



The distribution of radiative forcings in these two sub- 
sets is shown relative to the cloud-free radiative forcing of 
41.3 Wm- 2 (Fig. [4}. The maximum radiative forcing from 
subset 1 is 40.2 Wm~ 2 ; all physically plausible cloud sets 
give a smaller radiative forcing than a cloud-free model. Sub- 
set 2 - of cloud sets which give Earth-like climate - has 
a mean radiative forcing of 34.6 with a standard deviation of 
1.3 Wm -2 . The radiative forcing from the cloud-free case 
is 4.9 standard deviations above the mean radiative forcing 
from realistic, Earth-like, clouds. 

Note that we perform these runs at the standard solar con- 
stant, as for all model runs herein. Given that CO2 is not a 
strong absorber in the shortwave spectral region, selecting a 
lower solar constant for the sensitivity experiment would not 



x 10 




20 25 30 35 
Radiative forcing (W m~ 



45 



Fig. 4. Histogram of radiative forcings from two subsets of cloud 
profiles. Cloud sets which give energy balance at the TOA (sub- 
set 1) in light grey, cloud sets which give energy balance at the 
TOA and are close to observed longwave and shortwave fluxes at 
the TOA superimposed darker (green online). Dashed vertical line 
(red online) shows radiative forcing cloud-free case for comparison. 



Table 3. Cloud properties used in case study. f u 



=0.66. 



Property 


High 


Mid 


Low 


Cloud top (hPa) 


250 


500 


700 


Cloud base (hPa) 


300 


600 


850 


Cloud fraction 


0.25 


0.25 


0.40 


Water path (g m -2 ) 


20 


25 


40 


Liquid or ice 


Ice 


Liquid 


Liquid 


Generalised effective size ( p.m) 


75 






Effective radius (p.m) 




1 1 


11 



cause any noticeable change to Fig. |4] For example, using a 
solar flux 80% of the present value yields forcings different 
by 0.2 to 0.3 Wm" 2 . 

3.4 Case study selection 

As discussed, there are many problems associated with se- 
lecting a set of cloud profiles. However, the radiative forc- 
ings from CO2 enhancement in all Earth-like cloud sets are 
closely grouped (Fig. |4]l and the mean of these is signifi- 
cantly different from the cloud-free case. This justifies defi- 
nition of a case study which can be used to represent Earth's 
clouds. To do this from subset 2, we additionally constrain 
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Fig. 5. Comparison of global annual mean energy budgets: two estimates for contemporary climate, (a) Trenberth et al. (20091, based on 
a composite of data and (b) Zhang et al. ( 2004), from ISCCP-FD data, compared to models used in this paper, (c) case study with real clouds 
(d) cloud-free model. 



cloud fractions (each layer and the resultant total) and water 
paths of each layer to be close to climatological values, op- 
timising for agreement with longwave and shortwave fluxes 
at the FOA. We found that, whilst shortwave fluxes could 
be found that were in good agreement with climatological 
values, the outgoing longwave fluxes were slightly too high 
in all cases from ensemble 2. Increasing the height of the 
clouds by 50hPa gives a better fit for longwave fluxes. Case 
study cloud properties are given in Table [3] and the radiative 
outcome in Fig. [5};. 

3.5 Cloud-free case 

In order to compare calculated cloudy and cloud-free radia- 
tive forcings we need a cloud-free model as a comparison 



case. To generate this, we follow [Kasting et aL ( 1984 1 and 
tune the surface albedo of the GAM profile to achieve energy 
balance at the top of the atmosphere for a clear sky profile. 
The required surface albedo is 0.264. 



4 Real clouds and cloud-free model compared 

First, consider the energy budget at standard conditions rel- 
ative to observational climatology (Fig. [5]). Our case study 
with real clouds (Fig. [5};) is in very good agreement with ob- 



servational climatology (Fig. |5^,b). By contrast, almost all 
of the variable fluxes in the cloud-free model (Fig. |5jl) are 
markedly different; omitting clouds means that the global 
energy budget is not properly represented. Overall in the 
cloud-free model, more absorption of solar radiation (only 
8 1 W m~ 2 of outgoing shortwave radiation is reflected rather 
than 106 Wm~ 2 , a lower overall planetary albedo) is bal- 
anced by a weaker greenhouse effect (with an elevated out- 
going longwave flux of 261 Wm~ 2 rather than 236 Wm -2 
and depressed downward longwave at the surface). 

Whilst no 1-D model can perfectly represent global cli- 
mate, our real cloud case study, which is constrained by ob- 
servational cloud climatology, gives good agreement with the 
observed energy budget. This justifies using it as an internal 
standard, against which the cloud-free model can be com- 
pared. 

Again at standard conditions, compare the spectrally re- 
solved fluxes (Fig.[6|l. In the shortwave, the difference in ad- 
sorption between cloud-free and real cloud models (Fig. [6^) 
has the same shape as the Planck function of solar radiation 
(Fig. [6}:). This is because the surface albedo is constant with 
wavelength by definition and the wavelength dependence of 
cloud scattering is weak. Rayleigh scattering is spectrally de- 
pendent (short wavelengths are preferentially scattered), but 
this is a small term (14.5 Wm -2 in the cloud-free case). By 
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Fig. 6. Comparison of spectrally resolved energy budgets in 
Real Cloud (RC, blue) and Cloud Free (CF, red) models. Black 
lines are for both cases. Green is for differences (CF— RC). 
(a) F iiw A f o r both cases in black, F t| w A in colours, (b) F tLw 
for both cases in black, F tSw A m colours. (c) Absorp- 
tion of solar radiation: A=F |lw A — F tsw A (d) Greenhouse 
effect: G=F tLw ~F tLw A (e) Difference in solar absorp- 
tion: D A =A(CF)-A(RC) (f) Difference in greenhouse effect: 
L> G =G(CF)-G(RC). 



contrast, in the longwave, there is strong spectral dependence 
in the differences between the real cloud and cloud-free mod- 
els. The cloud-free model has a weaker greenhouse effect 
than real clouds in the water vapour window region. Whilst 
other spectral regions are optically thick (with gaseous ab- 
sorption by water vapour and carbon dioxide dominating) the 
water vapour window is optically thin and the cloud green- 
house is important. 

Now consider the effect of changing CO2 concentration 
(Fig. |7J. Radiative forcing is strongly overestimated by the 
cloud-free model relative to our real cloud case study; to pro- 
duce a given radiative forcing, twice as much CO2 is needed 
with the real cloud case study than is indicated by the cloud- 
free model. 

The radiative forcing in the longwave is an order of magni- 
tude larger than the radiative forcing in the shortwave region, 
so we focus on the longwave region when comparing spec- 
trally resolved forcing (Fig. IS). The greenhouse effect with 
real clouds is stronger at standard conditions than the cloud- 
free model (inclusion of cloud greenhouse). However, the 
greenhouse forcing (increase in strength of the greenhouse 
effect), is larger in the cloud-free model. This is true across 
all bands where CO2 imparts a greenhouse effect and is most 
important in the water vapour window. Here, the atmosphere 




10 10" 
pC02 (bar) 



10" 



Fig. 7. Radiative forcing with increasing pC02. Real clouds in blue 
and cloud-free in red. Present pC02 marked (*). 




10 20 

Wavelength, X (fim) 



10 20 
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Fig. 8. Comparison of spectrally resolved longwave forcings for in- 
crease from standard to 50000ppmv CO2 in real cloud (RC, blue) 
and cloud-free (CF, red) models. Green is for differences in cases 
(CF— RC) and black is for fluxes common between cases, (a) RC: 
FtLw in black, Ftw* dashed blue for standard C0 2 and solid 
blue for elevated CO2; (b) CF: F tLw in black, F t£ w A dashed red 
for standard CO2 and solid red for elevated CO2; (c) greenhouse 
forcing from increased CO a : G=G(HighC0 2 )-G(StdC0 2 ) 
where G=F , tLw f — ^Tlw A ; (d) difference in greenhouse forcing: 
De=G(CF)-G(RC). 



10 



C. Goldblatt and K. J. Zahnle: Clouds and the Faint Young Sun Paradox 




5 10 15 20 25 30 

Wavelength (/im) 



Fig. 9. Absorption cross sections for CO2 (green) and water vapour (purple) from HITRAN (shown at 900 hPa and 285 K). Horizontal lines 
indicate the cross section for which the gas has an optical depth of unity, solid purple for the GAM water vapour column, dashed green for 
various CO2 concentrations. The column depth of the atmosphere is 2.1x 10 25 molecules cm -2 . 



is optically thin in the absence of clouds, so the effect of in- 
creasing CO2 is large even though its absorption lines are 
weak (Fig. |9j. With clouds, these regions will be optically 
thicker initially so increasing CO2 has less of an effect. 

At 15 [am, increasing CO2 causes increased longwave 
emission. This is due to increased emission in the strato- 
sphere and is therefore unaffected by tropospheric clouds. 

Our GAM profile includes O3 which absorbs at 9.5 um and 
9.7 um. This would be absent in the anoxic Archean atmo- 
sphere, making the water vapour window optically thinner. 
The overestimation of forcing by the cloud-free model is, 
therefore, likely larger than suggested here and even more 
CO2 would actually be needed to cause equivalent warming. 

The other perturbation to consider is the change in in- 
coming solar flux. The cloud-free model absorbs a higher 
proportion of the incoming solar flux (has a lower plane- 
tary albedo), so will have a proportionately larger response 
to changing solar flux. For an 20% decrease in solar flux, 
representative of the late Archean, the decrease in absorbed 
solar flux is 52.3 Wm~ 2 for the cloud-free model compared 
to 47.2 Wm -2 for the real-cloud case study. Taking the 
Archean to be lower solar flux but higher CO2, this error 
is of opposite sign to the error in radiative forcing from in- 
creased CO2 and around half the magnitude. Whilst these 



errors could be said to partially offset in these conditions, 
reliance on errors of opposing sign is not strong. 

5 Variation of cloud and surface properties 

The problem of cloud feedback on climate change is noto- 
riously difficult. We do not attempt to address this in full; 
rather, we explore how variations in cloud amounts and prop- 
erties could affect climate. In all cases here, our baseline case 
is the real cloud case study and we consider the radiative ef- 
fect of changes in cloud or surface properties. As comparison 
values, if we increase or decrease the humidity in the model 
profile by 10% (50%), the radiative forcings are 1.7 Wm -2 
(7.8 W m- 2 ) and - 1 .9 W trr 2 (- 1 1 .2 W m" 2 ), respectively. 
A fainter sun in the late Archean is equivalent to a forcing of 
around — 50Wm~ 2 (assuming a planetary albedo of 0.3). 

5.1 Surface albedo 

In the cloud-free model, the use of a non-physical surface 
albedo means that real changes in surface albedo cannot be 
considered. This limitation is removed with explicit clouds. 
We limit discussion here to changes in surface albedo not 
from ice, though the use of a physically realistic surface with 
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albedo 



Fig. 10. Radiative forcing with changed surface albedo. (*) is the 
case study and ( o ) is the end-member case of an ocean covered 
planet. 



RC means that a parameterised ice-albedo feedback could be 
included in 1-D climate models, a significant improvement 
on the status quo. 

The surface albedo we use of 0.125 represents a weighted 



average of land (0.214) and ocean (0.090) albedos (Trenberth 



et al. 2009). Continental volume is generally thought to have 
increased over time, with perhaps up to 5% of the present 
amount at the beginning of the Archean and 20%-60% of the 



present amount by the end of the Archean (e.g. |Hawke s worth 



and Kemp 2006). In Fig. 10 we consider a range of varia- 
tion of surface albedo appropriate for a changed land frac- 
tion. For the end-member case relevant to the Archean of 
a water- world, the radiative forcing is 4.8 Wm~ 2 . Without 
land, relative humidity would likely be higher, contributing 
extra forcing. 

5.2 Cloud fraction and water path 

There are more clouds over ocean than land (Fig. [2|. The 
zonally uninterrupted Southern Ocean is especially cloudy. 
One might therefore expect that when there was less land 
there would have been more cloud, and more still if there was 
a greater extent of zonally uninterrupted ocean. Comparison 
of the Northern and Southern Hemispheres of Earth (Fig. |3J, 
the former having a higher land fraction, may be indicative 
of the minimum expected degree of variation. The Southern 
Hemisphere has 20-50% greater cloud fraction in each layer 
than the Northern Hemisphere. 

We consider a wide range of water paths, from optically 
thin to thick clouds, and fractional cloud cover from zero to 



1 for each cloud layer. In Fig. 1 1 we show the radiative forc- 
ing from these clouds relative to no cloud in the given layer. 
The competing shortwave and longwave effects of changing 
clouds can readily be seen. Increasing fraction or water path 
causes a negative forcing in the shortwave region (more re- 



flection) but a positive forcing in the longwave (greenhouse 
effect). The greenhouse effect operates by absorption of ther- 
mal radiation emitted by a warm surface followed by emis- 
sion at a lower temperature. Therefore the magnitude of 
changes in the greenhouse effect varies with cloud height, 
as higher clouds are colder. For low and mid level clouds, 
shortwave effects dominate and increasing cloud fraction 
or thickness will cause a net negative forcing (cooling the 
planet). For high clouds, shortwave and longwave effects are 
of similar magnitudes so the character of the net response is 
more complicated. For water paths less than 350 g m~ 2 , high 
clouds cause a net positive forcing (greenhouse warming the 
planet). The converse is true above 350 gm -2 , but such high 
water paths would typically correspond to deep convective 
clouds, not high clouds (cirrus or cirrostratus) (|Rossow and| 



Schiffer 



1999[ ). Positive forcing is maximum for ~70 g m 



high clouds. 
5.3 Cloud particle size 

Cloud particle size depends very strongly on the availabil- 
ity of cloud condensation nuclei (CCN). Whilst the global 
mean droplet size is 1 1 um, this is biased by smaller droplets 
over land (average 8.5 um), where there are more CCN than 
over the ocean (average 12.5 pm). Over the ocean, around 
half of CCNs are presently derived from oxidation products 
of biogenic dimethyl sulphide (DMS), especially sulphuric 
acid (there are various oxidation pathways of DMS (e.g. 

2004), but only sulphuric acid can 



von Glasow and Crutzen 



cause nucleation of new droplets ( |Kreidenweis and Seinfeld 
1988}). The climatic feedbacks involving DMS ( |Charlson| 



et al.| 1987[ ) have been subject of long debate. Whilst DMS is 
prevalent today due to production by eukaryotes, other bio- 
genic sulphur gases are produced by bacteria, in particular 
hydrogen sulphide (H2S) and methyl mercaptan (CH3SH) 



(Kettle et al. 200Tb. These will react chemically to form sul- 



phates, which will provide CCN. 

We do not delve deeply into CCN feedbacks here, but 
accept that various changes in the Earth system (e.g. atmo- 
spheric oxidation state, sulphur cycle, volcanic fluxes, bi- 
ological fluxes) may well have changed CCN availability. 
Fewer CCN give larger cloud drops, which should both rain 
out quicker (so less cloud) and be less reflective. Conversely, 
more CCN give more extensive and more reflective clouds. 

We consider the effect of changing liquid droplet size 
by factors of 0.5, 1.5 and 2 relative to the case study 
(r e jf=ll um) and ice particle size by factors of 0.5, 1.5 and 
1.87 relative to the case study (Z?GE=75um; the maximum 
of the parameterisation used is 



12 



we 



140 um). In Fig 
show the net (shortwave plus longwave) radiative forcing 
from changing particle size for all water paths and fractions. 
The effect is strongest for low clouds. With no change to 
cloud fraction or water path, increasing r e ff by 50% gives 
a forcing of 7.5 Wm~ 2 and doubling r e g gives a forcing 
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Fig. 11. Cloud radiative forcing with cloud fraction and water path relative to no cloud in that layer. For each cloud layer, these properties are 
varied whilst the clouds in other layers remain fixed at the case study values, marked (■#-). Particle sizes, ice/water ratio and height are as case 
study. Colour/contour scale is in Wm~ 2 . For comparison, resolution of the late Archean FYSP would require a forcing of approximately 
50Wm -2 . 



of 10.4 Wm . Decreasing r e g by 50% gives a forcing of 
-13.6Wm- 2 . 

Satellite observations of the modern ocean ( |Breon et al.[ 
2002) suggests a limit on how large droplet size actually be- 
comes in nature. Particle size is rarely larger than 15 pm, 
even in the remotest and least productive regions of the 
ocean. Here, the DMS flux is low and remaining CCN derive 
from abiological sources (e.g. sea spray). r c ff=15iim can 
then be seen as the baseline case for lower CCN availability, 
corresponding to a 36% size increase relative to present day 
mean (20% relative to present day ocean). 

If there was a larger CCN flux, the droplet size for clouds 
over land (r e g =8.5 iim, 23% less than mean) is an indicator 
of likely droplet size. 

Larger droplets will rain out more effectively, but model 



representations of this feedback vary dramatically (Pen- 



|ner et al.| [20061 l K " m P and Pollard| |2008). For the case 
of r e ff=17|am droplets over the ocean, |Kump and Pollard] 
(2008) choose a mid-strength assumption of this feedback, 
implying a decrease of water path by a factor of 2.2. This 



is marked (x) in the low cloud, 16.5 pm panel of Fig. 12 
the radiative forcing is then 15.4 Wm~ 2 , twice that of solely 
increasing droplet size. Clearly, an increased precipitation 
feedback is of first order importance and must be treated 
carefully in any model addressing the climatic effect of 
changed particle size. 

5.4 Cloud height 

To test the sensitivity to cloud height, each cloud layer is 
raised or lowered lOOhPa and the forcing is calculated rel- 
ative to the standard heights for all water paths and frac- 
tions. As temperature decreases with height, higher clouds 
emit at a lower temperature. It is this longwave effect which 
is dominant. There are only small changes in the shortwave 
effect, due to changed path length above the cloud (a greater 
path length means decreased insolation due to more Rayleigh 
scattering in the overlying atmosphere). In Fig. [13] we show 
only the net forcing. For the high clouds in the case study, 
which cover one-quarter of the sky and have a water path 
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Fig. 12. Change in net cloud radiative forcing with cloud particle size, across a range of cloud fractions and water paths. Particle size is varied 
for each layer independently (values in subplot titles), whilst all properties for other cloud layers remain as case study. The change is shown 
relative to the case study (so panels for r c ff =11 (jm and Dge=75 um contain the same information as Fig.[TT|net fluxes). Cloud fraction and 
water path for case study are marked (*); values at the point of these markers are for changing particle size only, values elsewhere in each 
panel are for changing water path or fraction too. Markers ( x ) and (+) refer to reduction in water path by factors of 2.2 and 3.7, respectively, 
for compar ison to|Rosing et al.](|2010), as discussed in the text. Marker (*) corresponds to the relatively thick and maximum extent clouds 
invoked by Rondanelli and Lindzen 1 2010 1. Colour/contour scale is in W m~ 2 . 



of 20 gm~ 2 , the effect of raising them is relatively small 
(2.8 Wm~ 2 ). There is a larger forcing from raising clouds 
which are thicker or cover more of the sky initially; the 
greater the radiative longwave effect of the cloud at its stan- 
dard height, the greater the effect of changing it height would 
be. 

Changes in the temperature-pressure structure of the at- 
mosphere might have induced changes in clouds. The 
Archean atmosphere was anoxic and did not have an ozone 



layer (Kasting and Donahue 1980 Goldblatt et al. 2006 1. 



Consequently, there would likely not have been a strong 
stratospheric temperature inversion, and deep atmospheric 
convection may have reached higher altitudes, where the at- 
mosphere is colder. A major source of high clouds is detrain- 
ment of cirrus from deep convective clouds. Where detrain- 
ment is due to wind sheer, this could then result in higher 
clouds. Conversely, without an inversion a the tropopause, 



cumulonimbus incus (anvil shaped clouds) will not form. As 
the forcing from raising the high clouds in the case study is 
small, other climatic effects might be larger (loss of ozone as 
a greenhouse effect and lower stratospheric emission temper- 
ature). Also, the pressure of Archean atmosphere was likely 
not 1 bar. Not only was there no oxygen (0.21 bar today), but 



the nitrogen inventory was likely different (Goldblatt et al. 



2009a I. Varying pressure would have changed both the lapse 



rate and tropopause pressure (Goldblatt et al. 2009a) 



6 Evaluating cloud-based proposals to resolve the Faint 
Young Sun Paradox 

6.1 Increased cirrus 



Rondanelli and Lindzen (20101 proposed that near total cov 
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Fig. 13. Cloud radiative forcing with changed cloud height relative 
to standard height clouds (see Fig. [TT| >. Colour/contour scale is in 
Wm~ 2 . 



erage of cirrus clouds could resolve the FYSP. Their pro- 
posed mechanism is that the planet would be colder and have 
lower sea surface temperatures, which would give more cir- 



rus coverage (the controversial "iris" hypothesis of Lindzen 
et al. , 2001), acting as a strong negative feedback on tem- 
perature. The first premise here, of colder temperatures, 
is contrary to the geological record; this suggests less fre- 
quent glaciation through the Archean and Proterozoic than 
in the Phanerozoic, not more. The second premise, of 
strong cloud feedback, is based on a statistical relationship 



for Earth's tropics ( jLindzen et al. |2001|l the authenticity of 



lartmann and Michelsen 



which has been questioned (e.g 

2002[ l. Application to very cold tern 



2002 Chambers et al. 



peratures requires an extreme and unverifiable extrapolation. 



Rondanelli and Lindzen ( 20 1 ) describe the high level clouds 
they use as "thin cirrus"; we note that the clouds they use ac- 
tually have twice the water path of our standard high clouds. 
In their sensitivity tests, using a thinner high level clouds 
gives a weaker effect. 

Here, we consider what would be required of cirrus or 
other high level clouds for them to resolve the FYSP. In- 
formed by the experiments above, we construct an optimum 
cirrus cloud for warming: relative to our case study we make 
it 3.5 times thicker (a water path of 70 gm~ 2 ) and make it 



gestion of Rondanelli and Lindzen 2010 1. This gives a forc- 
ing of 29.0 Wm -2 , insufficient to counter the ~50Wm" 2 
deficit from the FYSP. If, in addition, we raise the cloud by 
lOOhPa (base at 200 hPa, making the cloud 14 K colder) the 
total radiative forcing becomes 50.7 W m~ 2 . 

In principle, high clouds can resolve the FYSP. In prac- 
tice, the requirement for total high level cloud cover seems 
implausible and the requirement that the clouds are higher 
(colder) is difficult to justify. That it takes an extreme end- 
member case to provide only just enough forcing to resolve 
the FYSP suggests that resolution with enhanced cirrus only 
is not a strong hypothesis. 

6.2 Decreased stratus 



|Rosing et al.|(j2010| propose that there were less CCN avail- 
able in the Archean, due to lower DMS emissions prior to the 
oxygenation of the atmosphere and widespread occurrence of 
eukarya. They suggest an increase in droplet size from 12 um 
to 20 or 30 pm. Even over unproductive regions of today's 
oceans, the effective radius of cloud particles rarely exceeds 
15 um (Breon et al. 2002| l, so it is difficult to see how such 
large effective radii could be justified. Larger droplets lead 
to more rain, so should make clouds thinner. To account for 



this, |Rosing et al. (2010) arbitrarily decrease the liquid wa- 
ter path of their stratus clouds by a factor of 3.7, which is at 
the high end of likely decreases (jPenner et al. 2006 1. Even 
with these very strong assumptions, their model temperature 
is continually below the present temperature before 2 Ga. 

In our framework of radiative forcings, the effects of 
changing effective radius and cloud water path are shown in 



Fig. 12 For the strong but arguably plausible case (discussed 
in Sect. |5.3[ l of doubling the effective radius and decreas- 
ing water path by a factor of 2.2 gives a radiative forcing of 
15.4 Wm~ 2 . For the yet stronger case of doubling the effec- 
tive radius from 1 1 um to 22 um and decreasing cloud water 
path by a factor of 3.7, the radiative forcing is 20.5 W mT 2 . 
Removing low cloud entirely gives a forcing of 25.3 W m~ 2 . 
We therefore conclude that reducing stratus cannot by itself 
resolve the FYSP. 

A separate hypothesis (Shaviv 2003 , Svens mark| |2007| l 
proposes less stratus on early Earth due to fewer galactic 
cosmic rays being incident on the lower troposphere. The 
underlying hypothesis is of a correlation between galactic 
cosmic ray incidence and stratus amount, through CCN cre- 
ation due to tropospheric ionization ( jSvensmark and Friis-| 
Christensen 1997 Svensmark"] |2007| l. This hypothesis has 
been refuted (e.g. Sun and Bradley| |2002| [Lockwoo d and 
Frohlich 2007, Kristjansson et al. 2008; Bailer-Jones 



|Calogovic et al.| |2010[ |Kulmala et al.| |2010 and references 



2009 



therein): galactic cosmic rays cause the formation of at most 
10% of CCNs and there is no correlation between galactic 
cosmic ray incidence and cloudiness. Also of note is that 



Shaviv (2003 1 requires a highly non-standard climate sen- 
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sitivity to force his model. In the sensitivity test where he 
uses a more standard climate sensitivity, it results in a mean 
temperature ~ 0°C during the Archean. Even if the underly- 
ing hypothesis had not been refuted, the same arguments as 
above would apply: plausible decreases in stratus are insuf- 
ficient to resolve the FYSP. 



7 Conclusions 

When calculating radiative forcing from increased green- 
house gas concentrations, we find that omitting clouds leads 
to a systematic overestimate relative to models in which 
clouds are included in a physically based manner. With 
0.1 bar CO2 (the relevant quantity for a CO2 based resolu- 
tion to the Faint Young Sun Paradox in the late Archean) 
the overestimate of radiative forcing from modelling with- 
out clouds approaches 10 Wm~ 2 in our model, equivalent 
to the clear-sky forcing from lOOppm CH4. As the radia- 
tive transfer code we use underestimates forcing from CO2 at 
this level, and we include O3 in our profile, the difference be- 
tween the real cloud case study and the cloud-free model here 
must be seen as a lower bound on the error from omitting 
clouds. For other greenhouse gases, especially those which 
absorb strongly in the water vapour window, the overestima- 
tion by a cloud-free model would likely be larger. This would 
affect calculations of the warming by methane and ammonia 
and of recently proposed Archean greenhouse gases, ethane 
( |Haqq-Misra et~aL]|2008l > and OCS ( |Ueno etal.||200^ . 

The question of what direct effect clouds might have is 
a more interesting and difficult one. We can address this best 
by considering what radiative forcing can be generated in 
both the shortwave and longwave spectral regions by chang- 
ing cloud physical properties, and whether such changes in 
cloud physical properties can be justified. 

For solar radiation (shortwave), low level stratus clouds 
have the greatest effect. Removing them from the model en- 
tirely gives a forcing of 25 Wm~ 2 . Even this end-member 
falls short of the 50 Wm~ 2 that is needed to resolve the 
FYSP. A more plausible combination of reduced fraction and 
water path and increased droplet size would give a maximum 
forcing of 10-15 Wm -2 . However, suitable justification for 
these changes does not come easily. jRosing et al. ( 20 1 0) > 
asserted that DMS fluxes would be low in the Archean, 
but there may well have been other biological and chemi- 
cal sources for the sulphuric acid on which water condenses 
(DMS is a precursor to this). For example, methyl mercap- 

|200T) . 



tan is produced abundantly by bacteria (Kettle et al. 



Observations of clouds show that the effective radius rarely 
becomes larger that 15 um (Breon et al. 20021, which im- 



plies that regionally low CCN flux does not lead to very large 
droplets. If there were less land early in Earth's history and 
more zonally uninterrupted ocean, one might expect there to 
be more cloud rather than less (similar to how there is greater 
cloud fraction in the Southern Hemisphere than the North- 



ern Hemisphere). Also relating to low cloud, Shaviv ( 2003| > 
and Svensmark ( 2007 1 contend that less galactic cosmic rays 
were incident on the troposphere during the Archean and this 
would have led to less stratus. However, the underlying hy- 
pothesis for this has been refuted. 

For terrestrial radiation (longwave), high level clouds are 
most important as they are coldest (the greenhouse effect de- 
pends on the temperature difference between the surface and 
the cloud). The end member case is 100% coverage of high 
clouds which are optimised for their greenhouse effect, be- 
ing both thicker and higher than our case study. Such an end 
member case gives a forcing of 50 Wm~ 2 , which would just 
be sufficient to resolve the FYSP. However, physical justifi- 
cation for any of the required changes is lacking. Ron danefE| 
|and Li ndzen (2010) invoke a controversial negative feedback 
of increased cirrus fraction with decreased temperature (the 



"iris" hypothesis of Lindzen et al. 2001 1, but a true resolu 



tion to the FYSP should give temperatures equal or higher 
than present. Thus, even if the "iris" hypothesis was cor- 
rect, it would act to oppose warming. It is difficult to think 
of other mechanisms to make high clouds wider and thicker. 
Whether clouds should have been higher in the Archean may 
warrant more study. The absence of the strong stratospheric 
temperature inversion presently caused by ozone might con- 
tribute. However, without increase in fraction or cloud water 
path, the forcing will likely be less than 5 W m~ 2 . 

The question then naturally arises: How should one model 
early Earth climate? Some would look first towards a gen- 
eral circulation model (GCM), in order to better represent the 
dynamics on which clouds depend. We disagree. Whilst dy- 
namics are certainly important, it is unrealistic to think that in 
the near future clouds could be resolved in a global scale cli- 
mate model applicable to palaeoclimate. Even in "high reso- 
lution" models used for anthropogenic global change, cloud 
processes are parameterised sub-grid scale. As one moves to- 
wards deep palaeoclimate research, one moves further from 
the present atmospheric state for which the model may have 
been designed and can be validated. A larger model therefore 
introduces greater, and harder to track, uncertainty. Consid- 
ering what radiative forcing or warming a given mixture of 
greenhouse gases will impart is a first-order question, and 
one which should be answerable with a first-order model. 
A 1-D model is sufficient for this, but clouds must be in- 
cluded. The appropriate starting point would likely be a 
model with fixed cloud optical depth and fixed cloud top 



temperature (see, for example |Reck 1979 Wang and Stone 



1980). For any proposed change to clouds, very great atten- 
tion is needed to the feasibility of the mechanism involved. 
To model these, one should probably look towards a cloud 
microphysics resolving model, coupled to appropriate mod- 
els of CCN supply and chemistry. 

A stronger greenhouse effect likely contributes the largest 
part of the forcing required to keep early Earth warm. It 
is important to remember, however, that the forcing from a 
greenhouse gas depends on the logarithm of its abundance. 
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Thus, a modest forcing from clouds could have a large effect 
on how much of a greenhouse gas in needed; indicatively, 
a 10 Wm~ 2 from CO2 requires a increase in concentration 
of a factor of 2 to 3. A large atmospheric CO2 reservoir 
(~ 0.1 bar) may be slow to accumulate, as geochemical pro- 
cesses (principally volcanic outgassing) contributing linearly 
to concentration, so any other forcings may be rather useful 
in resolving the FYSP. 

In summary, it is necessary to include clouds in climate 
models if these are to be accurate. Resolution of the faint 
young sun paradox likely requires a combination of a few 
different warming mechanisms, including strong contribu- 
tions from one or more greenhouse gases. Changed clouds 
could contribute warming, but this has yet to be justified - 
and cooling caused by cloud changes is equally possible. Fu- 
ture work will no doubt propose novel mechanisms to change 
clouds. We hope that the results presented here will facil- 
itate quick and accurate look-up the climatic effect of such 
changes. Proposed cloud-based resolutions with only limited 
greenhouse enhancement are not plausible. 
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